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ABSTRACT 


The  effects  of  the  vibrational  excitation  and  deexcitation  of  nitrogen 
molecules  in  an  electrical  discharge  on  the  electron  energy  distributions  and 
transport  coefficients  are  investigated  theoretically.  The  electrons  collide 
with  vibrationally  excited  molecules  either  gain  or  lose  energy  which  results 
in  the  redistribution  of  the  electron  energy.  The  distribution  function  is 
calculated  numerically  by  solving  the  Boltzmann  equation,  using  a  complete  set 
of  elastic  momentum  transfer  and  inelastic  cross  sections.  The  presence  of 
vibrationally  excited  molecules  which  collide  and  exchange  energy  with  the 
electrons  has  a  strong  effect  on  thq  electron  energy  distribution  function. 
This  effect  is  demonstrated  for  various  vibrational  temperatures.  The  energy 
distributions  are  calculated  at  the  ratio  of  electric  field  to  gas  density  E/N 
which  ranges  from  10"iD  to  20  x  10  V  cm  .  Generally,  as  the  vibrational 
temperature  increases,  the  electrons  become  richer  in  the  high-energy  tail  of 
the  distribution  function. 


I.  INTRODUCTION 


The  evolution  of  the  electrical  conductivity  generated  by  a  propagating 
relativistic  charged  particle  beam  in  a  neutral  gaseous  medium  has  been  an 
interesting  subject  for  many  years.  The  ionization  of  the  gas  by  charged 
particle  beams  generates  secondary  electrons  which  either  lose  or  gain  the 
energy  through  elastic  and  inelastic  collisions  with  the  gas  species.  The 
conductivity  of  the  plasma  channel  can  be  determined  if  one  knows  the  electron 
energy  distribution,  which  results  from  the  interaction  of  the  plasma  with 
molecules  and  beam  induced  electromagnetic  fields.  For  simplicity,  it  is 
often  convenient  to  assume  that  the  electron  velocity  distribution  is 
Maxwellian.  However,  the  degree  of  accuracy  of  the  Maxwellian  assumption 
depends  on  the  plasma  conditions  and  the  magnitude  of  E/N. 

Numerical  calculations  of  the  electron  energy  distribution  in  typical 
electrical  discharge  in  Ng  have  been  studied  by  many  authors1”5  by  solving  the 
Boltzmann  equation.  Nitrogen  was  chosen  because  of  its  high  abundance  in  the 
atmosphere.  Besides,  the  cross  sections  for  the  elementary  processes  in 
nitrogen  have  been  studied  and  well  documented.  Once  the  energy  distribution 
function  is  obtained,  the  transport  coefficients,  such  as  the  electron  drift 
velocity,  the  ratio  of  the  diffusion  coefficient  to  the  electron  mobility  D/u 
(the  characteristic  energy)  etc.  are  obtained  in  the  usual  manner.  In  general 
the  results  show  that  the  calculated  transport  coefficients  are  in  good 
agreement  with  the  experiment.^"^  Most  of  all,  the  calculated  distribution 

1  C 

functions  are  found  to  be  markedly  non-Maxwell i an. 

In  most  practical  cases,  however,  the  fractional  population  of  the 
vibrationally  excited  molecules  in  an  electric  discharge  can  be  very 
significant,  while  the  density  of  electronically  excited  molecules  is  still 

g 

low.  Electron  collisions  with  the  vibrationally  excited  molelcules  can 
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result  in  the  excitation  (or  deexcitation)  to  higher  (or  lower)  vibrational 
levels.  As  the  vibrational  temperature  increases,  the  molecules  have  a 
substantial  degree  of  vibrational  excitation.  It  is  the  purpose  of  this  paper 
to  take  into  account  the  presence  of  highly  populated  vibrational ly  excited 
molecules  which  modify  the  electron  energy  distribution  because  of  the 
electron  inelastic  and  superelastic  collisions  with  the  molecules.  The  aim  of 
this  paper  is  to  demonstrate  that  one  must  account  for  superelastic  collisions 
with  the  vibrational ly  excited  molecules  to  obtain  the  proper  electron  energy 
distribution.  In  certain  cases  one  must  also  consider  the  superelastic 
collisions  with  the  electronic  states.  However,  this  will  be  the  subject  of  a 
future  paper. 


3 


II.  THEORY 


In  the  electron  energy  distribution  calculationss  the  conventional  two- 
term  spherical  harmonic  expansion1  is  chosen  to  represent  the  angular 
dependence  of  the  electron  distribution  in  the  velocity  space.  Accordingly 
the  Boltzmann  equation  for  a  spatially  independent  electron  swarm  in  a 
spatially  uniform  gas  under  the  influence  of  a  steady  electric  field  E  can  be 
expressed  as1 


d  r  E2u  dfptu)1  2m  d  r  2n  t  *  KT  df<^n 

rin  1 .  ?  "  ”  du  M  du  u  ^m  ^  0  e  'j 


dU  L3N2Qm(u) 


M  du 


du 


(i) 


=  [(u  +  u*k1  f° (u  +  u*kl  x*°*k(u  +  u*k) '  uf°(u)  xAk(u)l 

■I  t(u '  uik’  -  V  XA* (u '  u«k}  ■  uto(u)  xAi(u,l 


where  fg(u)  is  the  isotropic  component  of  the  electron  energy  distribution 
function  with  electron  energy  u  expressed  in  electron  volts,  N  the  gas 
density,  K  the  Boltzman  constant,  e  the  electron  charge  and  T  the  gas 
temperature.  M  and  m  are  the  mass  of  the  gas  molecule  and  the  electron, 
respectively.  Qm  is  the  electron  -  Ng  momentum  transfer  cross  section. 

0Ak  refers  to  vibrational  and  electronic  excitation  and  ionization  cross 
sections,  u^k  is  the  characteristic  energy  loss  for  each  process,  and  is 
the  fractional  concentration  of  the  Jtth  species.  The  second  term  on  the  right 
hand  side  will  be  retained  only  if  the  corresponding  deexcitation  is  present 
and  is  the  cross  section.  Equation  (1)  can  be  reduced  by  integrating  once 
with  respect  to  u  which  yields  the  result 


4 


(2) 
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- — +  ~u  Vu>  Cf0^)  +  e  du  ^ 


3IT  om(u) 


u+u, 


'Ak  u 

m  “l  I  ufOXAQAkdu  +  ^  /  ufOXAQkAdu 
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where  u^  is  defined  as  the  energy  loss  for  a  collision  in  which  the  molecule 
changes  its  vibrational  state  from  A  to  k.  Equation  (2)  is  a  first-order 
integrodifferential  equation  which  can  be  solved  numerically.  The  numerical 
method  described  in  Appendix  V  of  Reference  1  is  used  for  the  present 
calculation. 


III.  NUMERICAL  RESULTS 

Let  us  consider  the  case  where  the  fractional  concentration  of 
vibrationally  excited  molecules  is  so  low  that  the  effect  of  their  presence  on 
the  electron  energy  distribution  can  be  neglected.  Consequently,  the  second 
term  on  the  right  hand  side  of  Eq.  (2)  is  omitted  because  no  deexcitation  is 
being  considered.  In  other  words,  the  nitrogen  molecules  are  all  in  the 
ground  state  so  that  =  1.  In  practice,  the  subscripts  A  can  be  eliminated 
entirely  in  Eq.  (2).  This  Boltzmann  transport  equation  can  be  integrated 
numerically  using  the  backward  prolongation  technique  of  Sherman.10 

The  processes  considered  include  rotational,  vibrational  and  electronic 

excitation,  ionization  and  dissociation.  For  vibrational  excitations  only  the 

first  eight  levels  are  considered.  The  energy  loss  by  the  ionization  and 

dissociate  ionization  is  treated  as  one  of  the  excitation  processes.  A 

complete  set  of  electron  cross  section  data  was  used  in  solving  the  Boltzmann 

equation  numerically  for  the  steady  state  condition.  The  electron  energy 

distributions  for  various  values  of  E/N  are  shown  in  Fig.  1.  The  E/N  values 

—16  7 

are  varied  from  1  -  20  x  10  V  cm  ,  corresponding  to  mean  electron  energies 
approximately  from  1  to  5  eV.  The  calculated  distributions  show  the  important 
features  of  the  electron  energy  loss  in  Ng.  Because  of  the  peak  value  of  the 
vibrational  cross  sections  near  2  eV,  it  is  difficult  for  the  electrons  to 
obtain  larger  energies  to  pass  over  the  vibrational  barrier.  Note  that  the 
distribution  function  fg(u)  has  been  normalized  as  in  Eq.  (3) 


For  each  distribution  function  the  calculated  drift  velocity  W,  the 
characteristic  energy  D/u  and  the  average  electron  energy  Eav  are  listed 
below: 


E/N  (V  cm2) 

1.0  x  10'16 
2.0  x  10-16 
5.0  x  10‘16 
1.0  x  10"15 
1.5  x  10”15 
2.0  x  10‘15 

If  the  Maxwellian 
the  energy  distributio 


W (cm/sec) 

D/pt 

1.79  x 

106 

0.98 

3.00  x 

106 

1.17 

6.27  x 

106 

1.38 

1.07  x 

107 

2.02 

1.48  x 

107 

2.96 

f— * 

• 

00 

X 

107 

3.74 

function  is  normalized  as 
becomes 


Eav^V> 

0.96 

1.04 

1.15 

2.23 

3.80 

5.02 

ndicated  by  Ea.  (3),  then 


f(u)  =  2% 


-1/2  T  -3/2 


-u/T 


where  Tg  is  the  average  electron  temperature  expressed  in  electron  volts. 
Therefore  the  Maxwellian  distribution  functions  for  any  Tg  would  appear  in  the 
log-linear  plot  as  a  straight  line.  For  comparison  we  plot  in  Fig.  2  the 
Maxwellian  functions  corresponding  to  each  curve  in  Fig.  1.  It  is  obvious 
that  the  calculated  electron  energy  distribution  functions  are  non- 
Maxwellian.  However,  it  seems  that  fg(u)  tend  to  be  more  Maxwellian  as  E/N 
increases.  Clearly  the  resulting  energy  distributions  have  shown  the 
important  features  of  the  electron -molecule  energy  exchange  processes  in  the 
electrically  excited  N2.  In  Fig.  3  we  show  the  calculated  electron  excitation 
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rate  coefficients  for  eight  ground  state  vibrational  levels  of  N2.  At  low 
E/N,  most  of  the  energy  is  dissipated  into  the  vibrational  excitations.  Also 
shown  in  Fig.  3  are  the  equivalent  excitation  rates  obtained  by  assuming  a 
Maxwellian  energy  distribution.11  It  can  be  seen  obviously  that  the 
Maxwellian  distribution  is  not  a  good  approximation  for  calculating  these  rate 
coefficients  until  E/N  exceeds  3  x  10“^  V  cm^. 

We  now  consider  the  influence  of  the  vibrationally  excited  molecules  on 
the  electron  energy  distribution.  If  the  molecules  have  a  substantial  degree 
of  vibrational  excitation,  one  must  take  into  account  the  superelastic 
collisions  of  the  electrons  with  the  vibrationally  excited  molecules.  In 
other  words,  the  excitation  and  deexcitation  should  be  considered  at  the  same 
time.  For  simplicity,  we  assume  that  the  vibrational  states  are  in  a 
Boltzmann  distribution  where  their  fractional  concentration  for  each 
vibrational  state  l  can  be  expressed  as 


vN 


total 


-0.3/T  -0.3I/T 

e  v)  e 


(5) 


where  Ty  is  the  vibrational  temperature.  Equation  (2)  has  to  be  solved  in 
full  with  the  right  hand  side  representing  collisions  of  the  first  and  second 
kind.  If  Ty  =  0  then  =  0  for  every  X  except  S.  =  0.  That  is  exactly  the 
case  of  no  vibrational  excitation  just  discussed.  The  entire  matrix  of  cross 
sections  for  vibrational  excitation  and  deexcitation  are  related  by  Eqs.  (6) 
and  (7) 


(6) 


Q0,i(u  +  “th1  +  “th1  ■  Qi,0(u)u  for  u  >  0 


(7) 
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where  0^  <  is  defined  as  the  cross  section  for  a  collision  in  which  the 

'  *J 

mulecule  changes  its  vibrational  state  from  i  to  j.  In  Eq.  (7),  use  has  been 
made  of  detailed  balancing  between  collisions  of  the  first  and  second  kind. 

is  the  threshold  energy  for  vibrational  excitation  of  the  nitrogen 
molecule  by  electron  impact  from  the  ground  state  to  ith  vibrational  state. 
For  calculations  of  the  distribution  function  with  vibrational ly  excited 
molecules,  the  contribution  to  the  right  hand  side  of  Eq.  (2)  due  to  the 
presence  of  the  vibrational  states  from  X  =  0  to  8  can  be  expressed  as 


Ik 


8  8 

•I  C-I  / 

A=0  k=A+l  u 


u+u 


xk 


ufoxAk 


Jt-1 

du  +  y 

k=0 


u-u 


ufoxA* 

xk 


du] 


where  the  two  terms  represent  the  excitation  and  deexcitation  of  vibrationally 
excited  molecules,  respectively.  Making  use  of  Eqs.  (6)  and  (7),  Eq.  (8) 
yields 


* 
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l 

k=J!+l 


U+UJlk  X-l  U+U£k 

/u  uf^du  +  Ju  uf0(u-uAk)QkJldu] 


(9) 


where  the  second  term  involves  the  distribution  function  fQ(u  -  uAk)  which  is 
unknown  and  causes  some  difficulty  to  use  the  backward  prolongation  technique 
to  solve  Eq.  (2).  Another  method  involving  a  matrix  inversion  is  necessary  to 
take  care  of  this  problem.  However,  since  the  second  term  can  be  omitted  at 
Ty  =  0,  we  can  at  least  predict  the  solution  of  fg(u)  at  Ty  =  0  using  the 
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background  prolongation  technique.  As  Tv  increases  graduallys  a  predictor- 
corrector  iterated  procedure  is  used  until  fg(u)  converges  to  steady  state 
solution. 

The  electron  energy  distribution  functions  are  calculated  for  various 
vibrational  temperatures  at  the  ratio  of  E/N  range  from  10"AD  to 
20  x  10“^®  V  cm^.  The  distribution  functions  for  E/N  =  3  x  10“^  V  cm^  are 
shown  in  Fig.  4  for  various  vibration  temperatures.  The  effect  of  the 
vibrational  temperature  on  the  distribution  starts  to  become  important  as  Ty  > 
0.04  eV.  As  the  vibrational  temperature  increases*  the  vibrational ly  excited 
molecules  transfer  their  excitational  energy  to  electrons  in  collisions  of  the 
second  kind,  so  that  the  electron  energy  distribution  becomes  relatively  rich 
in  the  high  energy  portion  of  the  distribution  function.  As  a  second  example, 
we  choose  E/N  =  5  x  10_i  V  cm  and  do  the  similar  calculations.  The  results 
are  shown  in  Fig.  5.  The  threshold  of  vibrational  temperature  for  effecting 
the  distribution  curve  is  around  0.05  eV.  As  E/N  increases,  it  seems  that 
higher  vibrational  temperature  is  necessary  to  increase  the  high-energy 
electrons.  Finally,  in  the  case  of  E/N  =  10  x  10“*®  V  cm^  the  distribution 
curves  are  shown  in  Fig.  6  with  the  vibrational  temperature  threshold  near 
0.1  eV.  As  Ty  increases,  there  is  a  decrease  in  fg(u)  in  the  0  -  2  eV 
range.  The  electrons  gain  a  lot  of  energy  from  the  vibrational  deexcitation 
especially  for  such  high  vibrational  temperature.  The  slow  electrons  gain 
sufficient  energy  to  pass  the  barrier  and  become  faster  electrons.  This 
results  in  an  increase  in  the  high-energy  portion  of  the  distribution 
function.  However,  as  Ty  continues  to  increase  to  0.8  eV,  the  distribution 
curve  near  2  eV  tends  to  be  more  Maxwellian  as  clearly  indicated  in  Fig.  6. 
Further  increase  in  Ty  beyond  0.8  eV  makes  minimal  change  in  the  energy 
distribution. 
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Finally,  as  to  the  dependence  of  the  transport  coefficients  on  the 
vibrational  temperature,  we  found  that  the  calculated  electron  drift 
velocities  remain  almost  unchanged  as  the  vibrational  temperature  increases. 
However,  the  characteristic  energy  and  the  average  electron  energy  are  strong 
functions  of  the  vibrational  temperature,  especially  for  higher  value  of 
E/N.  As  an  example,  for  E/N  =  1(T*5  V  cm2  Fig.  7  shows  the  calculated 
characteristic  energy  D/u  and  the  average  electron  energy  Eay  as  a  function  of 
the  vibrational  temperature.  The  results  indicate  that  vibrational 
temperature  variations  have  an  important  influence  on  the  high-energy  portion 
of  the  electron  distribution  in  electrically  excited  gases  and  thus  on 
electronic  excitation  and  ionization.  Therefore  the  effect  of  the  vibrational 
excitation  and  deexcitation  of  molecules  on  the  electron  energy  distributions 
and  transport  coefficients  can  not  be  ignored. 
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The  electron  energy  distribution  functions  and  transport  coefficients  in 
molecular  nitrogen  are  calculated  for  the  moderate  value  of  E/N  range  from 
1  x  10"xo  to  2  x  10"1  V  cm.  The  calculations  are  carried  out  by  numerically 
solving  the  Boltzmann  equation,  using  a  set  of  accurate  electron-collision 
cross  sections.  The  results  indicate  that  the  electron  distribution  is  not 
Maxwellian  and  the  corresponding  electron  impact  excitation  rate  coefficients 
for  eight  ground  state  vibrational  levels  are  significantly  different  from 
Maxwellian  rates,  especially  for  low  E/N. 

The  emphasis  of  this  study  has  been  placed  on  the  calculation  of  the 
distribution  function  in  the  presence  of  vibrationally  excited  molecules  and 
the  dependence  of  the  electron  energy  distribution  on  the  variation  of 
vibrational  temperature.  The  excitation  as  well  as  deexcitation  processes 
have  to  be  considered  in  order  to  obtain  the  true  electron  distribution 
function  as  our  calculations  clearly  indicate. 
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The  electron  energy  distributions  for  various  values  of  E/N. 

The  Maxwellian  functions  corresponding  to  each  curve  in 
Fig.  1. 

The  calculated  electron  excitation  rate  coefficients  for 
eight  ground  state  vibrational  levels  of  N£  are  compared 
to  that  obtained  using  a  Maxwellian  electron  energy 
distribution. 

The  electron  energy  distribution  funtions  for  various 
temperatures  when  E/N  *  3  x  10“  ^  V  cm^. 

The  electron  energy  distribution  functions  for  various 

1C  p 

vibrational  temperatures  when  E/N  =  5  x  10"A  V  cm  . 

The  electron  energy  distribution  functions  for  various 

IE  p 

vibrational  temperatures  when  E/N  *  10  V  cm  . 

The  calculated  characteristic  energy  D/p  and  the  average 
electron  energy  Eav  as  a  function  of  vibrational  temperature 
when  E/N  -  10“^  V  cm^. 
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